# Both lambda and delta move

max_W = 3000
W_values = seq(-max_W,max_W,by=10)
decalage = 0

lambdas = seq(4,2,length=20)
lambdas  = 10**-lambdas

deltas = seq(4,2,length=20)
deltas  = 10**-deltas

for (j in 1:length(lambdas)){
  
      grass = sigmo_pluie(W_values,lambdas[j], decalage = decalage)
      plot(
           W_values,grass,
           xlim=c(-max_W,max_W),
            ylim = c(0,1),
           main=paste("j=",j,"lambda=",round(lambdas[j],5))
           )
      abline(v=0,lty=2)
      
      grass = sigmo_coeff_competition(W_values,deltas[10], decalage = decalage)
      plot(
           W_values,grass,
           xlim=c(-max_W,max_W),
           ylim = c(0,1),
           main=paste("i=10",round(deltas[10],5))
           )
      abline(v=0,lty=2)
      
      ##############################################################################################
      
      palette <- rainbow(length(deltas))

      par(mfrow = c(1, 2))
      plot.new()
      legend("center", legend = paste("i=",1:length(deltas),"delta=",round(lambdas,4)), 
             col = palette, lty = 1,xpd = TRUE)
      
      plot(
           x=c(0),
           y=c(0),
           xlim=c(-max_W,max_W),
           ylim = c(0,0.5),
           xlab = "W", ylab = "",
           main=paste("j=",j,"lambda=",round(lambdas[j],5))
           )
      
      for(i in 1:length(deltas)){
         grass = sigmo_pluie(W_values,lambdas[j],decalage = decalage)*sigmo_coeff_competition(W_values,deltas[i],decalage = decalage)
         lines(W_values, grass, col = palette[i])
         abline(v=0,lty=2)
         abline(h=1,lty=2)
         abline(h=0.25,lty=2)
      }
}